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ABSTRACT 

A single degree of freedom aeroelastic computer model, WMSTAB 3 , has been employed to perform a parametric 
analysis of^HAWT blade behavior during yaw maneuvers. Over 1,000 different combinations of 63 and normal 
^equency were analyzed. 

The effect of 63 and flapping stiffness on flapping frequency, phase, and magnitude are discussed. The 

moments transmitted to the fixed system during yaw maneuvers are calculated and reduced to time constants 
of response to step changes in wind direction. The significance of the time constants for the configu- 
rations considered relative to yaw response rate and lag angle is discussed, along with their possible 
significance for large HAWT. 


: INTRODUCTION 

The 63 hinge Is a device which mechanically cou- 
ples rotor blade pitching and flapping. The 6 ^ 

hinge j s’ typically employed to stabilize lifting 
rotors^ by “increasing the flapping stiffness and 
reducing flapping. The stiffness increase also 
in creases the flapping frequency, thereby allowing 
flexibility in blade tuning to avoid possible 
structural resonances] 


Negative 63 is destabilizing. The effect of neg- 
atT^ 1 J is to reduce the flapping frequency, by 

Virtue of reduction in flapping stiffness. The 
re duction in f lapping frequency will be accom- 
panied by increasing flapping as long as the 
flapping frequency Ts approaching 1.0 per rev. 

The phase of flapping response to yaw direction 
and rate excitations also changes with adjustment 
of frequency. 

These qualitative observations of the effect of 
03 on the magnitude and phase of flapping 

response to yaw inputs suggested a quantitative 
analysis to determine optimum operating configu- 
rations for spring restrained teetering type wind 
turbines. {The 63 hinge is easily incorporated 

in such a rotor by canting the teeter axis and 
the normal to the feathering axis, as In Fig- 
ure 1 .) Because of the large range of flapping 
frequencies and candidate 63 angles considered, a 

simple aeroelastic model, WMSTAB 3 , was developed 
which could be rapidly executed, yet which 
retained sufficient sophistication to be accep- 
table. More than 1,000 possible operating con- 
figurations were ultimately analyzed to develop 
the data reported herein. 


or phase, or yaw performance, quantified by 
response time and/or restoring or damping 
moments. The operating conditions required for 
optimization are often contradictory. For 
example, it is possible to eliminate flapping 
In the direction of the tower due to either yaw 
position or rate excitations, but not both sim- 
ultaneously. 


ANALYTIC BACKGROUND 


The equations of motion were derived from funda- 
mental mechanics and a linear aerodynamic repre- 
sentation of the forcing function, paralleling 
Stoddard (Reference 1 ). The equation of motion 
of the flapping wind turbine blade is 
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The analysis indicates that it is possible to 
optimize rotor flapping with respect to magnitude 
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*The work reported herein was done as part of a 
Kaman Aerospace Corporation in-house R & D 
program. 
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into the eg^uation of motion, neglecting products 
of 3> a^, Uq and their derivatives, and rewriting 

the equation of motion in the frequency domain, 
we have 


K + ^ Tan 5^ 


6 1 




IL Tan 6-3 


”2" - -y Uq Tan 5 ^ K” 1 ^3^®^ ^3 


A 




^^0 

JT 




nc 


ns 


1 , 1 '"o 

T ^0 ■ ■ T r* ' 


i 






where terms in harmonics of f greater than the 
fundamental are neglected. This system of equa- 
tions is solved in program WMSTAB. 

The use of linear aerodynamics for the forcing 
function retains the essential features of a real 
wind turbine blade, including pitch angle, built- 
in twist, induced velocity, and blade mass, while 
sacrificing the accuracy of a more comprehensive 
aerodynamics analysis. Similarly, the single 
degree of freedom, flapping, retains the essential 
feature of a more comprehensive aeroelastic code, 
the mobility of the blade, while minimizing the 
computational difficulty. Because of the compu- 
tational ease of the analysis, it is very useful 
for the reduction of a large number of candidate 
operating configurations into general trends. 

The program WMSTAB uses the structural character- 
istics of the Kaman 40 kW wind turbine, where 
applicable, as the baseline. Those variables 
dependent on rotor size are proportioned for 
other radii as follows: 


R « R^ 



e « 


Three rotor radii were analyzed to illustrate 
characteristics dependent on rotor size. These 
radii (4.88 m [16 ft]; 9.76 m [32 ft]; and 
19.51 m [64 ft]) were taken to be representative 
of small, intermediate, and large wind turbines, 
respectively. 

The equation of motion does not include an offset 
between the yaw axis and the rotor. This offset 
will have both aerodynamic and dynamic consequen- 
ces, principally arising from three sources. 

First, yaw motions impart an out-of-plane velocity 
component to the blade, given by 

^=axr 

where the underscore denotes vector quantities. 

For small values of yaw offset, the magnitude of £ 
is essentially equal to the blade radius at points 
where the aerodynamic forces are large. Hence, 
the effect of yaw offset on blade translations is 
second order, and is neglected. 

Second, for very low normal frequencies, the tilt 
of the thrust vector will dominate yaw dynamics. 
For a true teetered rotor, there are no other yaw 
restoring or damping moments than those given by 
the product of the in-plane component of rotor 
thrust and the yaw offset. The magnitude of these 
moments depends on the exact value of the yaw 
offset so they cannot be calculated for a general 
case. 

Third, the yaw offset will cause the inertia of 
the yawing system to be greater than the flapping 
inertia of the blades alone. The magnitude of the 
mass transfer effect again depends upon the magni- 
tude of the yaw offset, and cannot be computed in 
general , 


FLAPPING ANALYSIS 

Because of the obvious difficulty in referring to 
both the actual flapping frequency and the stiff- 
ness of the flapping spring as some multiple of 
as is customary, the term "normal frequency" was 
coined to refer to the flapping frequency of the 
blade when 6^ = 0. Thus, increasing normal fre- 
quency refers to increasing spring stiffness, and 
a given normal frequency blade may have either 
increasing or decreasing flapping frequency, 
depending on 6^. 

Flapping is defined as the amplitude of the har- 
monic portion of the B deflections. WMSTAB3 
includes only first harmonic flapping. It is 
probable that this harmonic would include the 
bulk of flapping energy in an actual wind turbine 
because the cyclic inputs leading to higher har- 
monic responses are generally lacking. 

Figure 2 shows the dependence of flapping fre- 
quency on normal frequency, which is determined 
by the stiffness of the flapping spring, and on 



6^ which greatly affects the strength of the 

aerodynamic spring. The range of normal frequen- 
cies was defined by practical considerations of 
rotor blade design. The relationship between 
decreasing flapping frequency and decreasing 6^, 

and the converse, is apparent. Flapping insta- 
bilities are suggested by the rapid decay of 
flapping frequencies in the neighborhood of large 
negative 6^ angles. 

Figures 3 and 4 show the dependence of total 
flapping induced by yaw position and rate, 
respectively, on normal frequency and ^ 

number of trends are evident. First, for <5^ ^ 0» 

increasing norm_al frequency decreases flapping. 
This Is net necessarily so for 6^ < 0. Second, 

the total flapping is determined more by the 
proximity of the flapping frequency to 1.0 P 
than by normal frequency. In fact, the amount of 
flapping available to any normal frequency, given 
suitable essentially constant. The third 

1t^ of major interest is that the flapping due to 
yaw position is independent of rotor size, but the 
flapping due to yaw rate increases linearly with 
rotor size for a fixed yaw rate. This is a conse- 
quence of the gyroscopic moments acting on the 
blade, and gives rise to the yaw dynamics effects 
discussed later. 

Figures 5 and 6 present the phase angle of 
flapping response to yaw position and rate, 
respectively. Horizontal blade position corres- 
ponds to ^ 90®. Consequently, flapping which 
takes place entirely in the horizontal plane has a 
phase angle of ^ 90®. Cewnparison of Figures 5 and 
6 indicates that it is not possible to totally 
eliminate flapping in the tower direction. Judi- 
cious selection of normal frequency and 6^ can, 

however, insure that flapping in the direction of 
the tower will always be small relative to total 
flapping. Note that these phase angles are inde- 
pendent of rotor size, and depend only on the 
flapping frequency (and the pitch angles, a second 
order effect). 

Finally, the flapping Induced by yaw direction 
will depend approximately linearly on wind speed, 
while the yaw rate induced flapping is essentially 
independent of wind speed. 

The data presented in this section may be used to 
estimate the yaw performance of a true teeter 
rotor (zero flapping spring) through the phase 
angles and flapping angles presented. Again, the 
results will depend on a particular design and are 
not included In this analysis, 

YAWING BEHAVIOR 

The yawing behavior of a wind turbine depends on 
three things; the aerodynamic spring rate, the 
total damping, and the yawing inertia. The gover- 
ning equation of motion 

^y’^y ^ Vy ^ moments 

can be solved as discussed 1n many texts of ele- 
mentary dynamics. The contribution of this study 


is the generation of linearized estimates for 
S ^y solution of the homogeneous 

equation of motion, after setting 1^ = Ij^, from 

which time constants of yaw response to step 
changes In wind direction were calculated. 

Figure 7 is an illustration of the aerodynamic 
spring constant, K^, determined by imposing a yaw 

angle of 10®, and normalizing the result to a yaw 
angle of 1.0 radian. The peak of the spring rate 
occurs at the 1.0 P flapping frequency for all of 
the normal frequencies shown. Instability is^ 
suggested by the very large spring constants in 
evidence for high normal frequencies and large 
negative values of 6^. The stabilizing effect of 

positive 6 2 is also evident by the convergence of 
spring constants for large positive 6^, sug- 
gesting that 6^ has come to dominate flapping 

dynamics. Note that the yaw spring constant is a 
linear function of rotor size. The spring rate is 
approximately linear with wind speed for fixed yaw 
angle. 

The damping constant is plotted in Figure 8. 

Much of what has just been described for Figure 7 
Is also true here. The most notable exceptions 
are that the damping rate was normalized to 
1.0 radi an/sec, and the peak damping occurs at 5^ 

angles smaller (more negative) than for the cor- 
responding spring rate. The yaw damping rate is 
Independent (to first order) of wind speed. 

The spring rates and damping rates are decoupled 
in this analysis. 

The yaw dynamics behavior can be summarized by 
calculating the time constant of response to a 
step change in wind (yaw) direction, as discussed 
above. These data are assembled in Figure 9. 

Two trends, increasing time constant first with 
Increasing 5^, and second, with increasing normal 

frequency, are immediately obvious. The time 
constants also increase linearly with rotor size. 

The time constant can be interpreted many ways. 
Also indicated on ordinate axes of Figure 9 are 
the time required to decay to 50^ amplitude, and 
the time required to decay to 5% amplitude, which 
numbers are linear transforms of the time con- 
stant. Yet another interpretation is that the 
product of the time constant and the rate of 
change of wind direction equals the lag angle 
between the wind direction and the axis of rota- 
tion. Because the time constants increase with 
Increasing rotor diameter, the data indicate that 
the lag angle for large rotors could be quite 
large, even If the rate of change of wind direc- 
tion is quite small. Since this lag angle 
degrades performance, the evaluation of system 
productivity, particularly for large rotors in a 
variable wind, should be evaluated with respect 
to the lag (yaw) angles almost certain to be 
present. 
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CONCLUDING REMARKS 

The incorporation of significant amounts of 6^ 

into a horizontal axis wind turbine blade greatly 
increases the options available to the designer. 
The 5^ hinge may be used to adjust magnitude or 

phase of the flapping response to either yaw rate 
or direction, thus controlling the mechanical and 
aerodynamic coupling with the fixed system 
degrees of freedom. The effect of 6^ on flapping 

frequency can also be exploited to optimize sys- 
tem dynamics. 

The reduction of the blade flapping data into 
yawing spring and damping constants allows the 
estimation of response time to changes in wind 
direction. The data indicate that large wind 
turbines may suffer large lag angles if operated 
where winds vary in direction. This may be a 
source of performance degradation previously 
overlooked. 
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NOMENCLATURE 

a lift slope = 5.73/radian 
C chord 

C^ yaw damping rate 

t cosine 
e off$et/R 

g gravitational acceleration 
flapping inertia 

yaw inertia 
flapping spring rate 
aerodynamic yaw spring rate 
Mj^ rotor blade mass 

P division by q 
q yaw rate 

q 

R rotor blade radius 

r distance from yaw axis to a point on 
the blade 

$ sine 

Uq crossw^nd/^^R 

out-of-plane velocity due to yaw rate 
Xg radius to c.g./R 


ocy yaw angle 

6 flapping angle 
3 q steady harmonic coefficient 

3^^ cosine harmonic coefficient 
3 -|^ sine harmonic coefficient 

4 

Y lock number ^ pCaR /l^ 

6^ angle between normal to feathering axis 
and teeter axis 

0Q built-in linear twist (total) 

0p pitch angle, + implies increasing thrust 
induced velocity/QR 
Uq wind velocity/aR 

p air density 
¥ azimuth angle 
n rotational speed 



Figure 1 - Teetered rotor showing 5^ due to 

offset between teeter axis and 
the normal to the feather axis. 
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QUESTIONS AND ANSWERS 


F,W,. Perkins 

Fromi 

Q; Can you physically explain the significance of 1.4 P, and its effect on yaw dynamics? 

A: 1,4 V is just that normal frequency which maximizes the product of Kq and which 

p 2 1> 

gives the yaw moment. Incerasing consequently the normal frequency y decreases 

0^^ more rapidly than increases . 

I 

I 
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